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ABSTRACT 

While in the standard cosmological model the accelerated expansion of the Universe is ex¬ 
plained by invoking the presence of the cosmological constant term, it is still unclear the true 
origin of this stunning observational fact. It is therefore interesting to explore alternatives to 
the simplest scenario, in particular by assuming a more general framework where the fluid 
responsible of the accelerated expansion is characterised by a time-dependant equation of 
state. Usually these models, dubbed dark energy models, are purely phenomenological, but 
in this work we concentrate on a theoretically justified model, the ghost dark energy model. 
Within the framework of the spherical collapse model, we evaluate effects of dark energy 
perturbations both at the linear and non-linear level and transfer these results into an observ¬ 
able quantity, the mass function, by speculatively taking into account contributions of dark 
energy to the mass of the halos. We showed that the growth rate is higher in ghost models and 
that perturbations enhance the number of structures with respect to the ACDM model, with 
stronger effects when the total mass takes into account dark energy clumps. 


1 INTRODUCTION 

In the standard cosmological model, the small initial fluctuations 
that were seeded during the phase of inflationary expansion, 
some 10“®® seconds after the Big Bang, are the origin of the 
large scale structures that we observe t o day, i.e., galax i es and 
cluste rs of galaxies dStarobinskvl Il98d : iH.GuthI Il98ll : iLindd 


of gravitational collanse Jcunn & Gotti 

I 972 I: Press & Schechter 

1974 IWhite & Reeslll978l: IPeeblesll 19931: ISheth & TormenI 

19991: 

Peacock! 1999: Barkana & LoebI 2001 

Peebles & Ratral 

2003; 

Ciardi & Eerraral I 2 OO 5 I: iBromm & Yoshidal l201lb. Most of the 


growth has taken place after the decoupling of photons and 
electrons. Th e sphe rical collapse model (SCM) introduced by 
iGunn & Gotj ( Il972h is a simple analytical model to study the 
evolution of the growth of cosmic structures. The scales of interest 
in the SCM are much smaller than the Hubble length and the 
velocities are non-relativistic. Therefore, pseudo-Newtonian grav¬ 
ity can be used to study the evolution of the overdensities. With 
pseudo-Newtonian gravity we refer to the fact that it is possible to 
use Newton’s hydrodynamical equations in an expanding Universe 
including relativistic contributions to the Poisson equation, i.e. the 
contribution of pressure terms. At early times, when the overden¬ 
sities are small, linear theory is able to follow the evolution of 
spherical overdense regions. In this phase, due to self gravity, the 
overdense region expands at a slower rate compared to the Hubble 
flow. At a particular scale factor, depending on the particular 
background history, the spherical region reaches a maximum 
radius and completely detaches from the background expansion. 
This is the so called tum-around epoch. Subsequently the collapse 
proceeds under the overdense region own gravity and the system 
reaches a final steady state with a specific radius due to virialization 


processes. The dynamics of the SCM depends strongly on the 
evolution of the background Hubble flow. The SCM was extended 
and improved in s everal works i Fillmgre_&Goldreich|_ 1^84 
Bertschingeil 198jj_ Hoffman & Shaha ml Il985l: iRvden & GunnI 
19871: lA vila-Ree se et al.l 1 19981: ISubramanian et al.l 2000l : 


Ascasibar et al. 20041: Williams et all 


20041) . More recently. 


this fo rmalism has been exten ded to include shear and rotation 
terms jdekPopolo et_^ |20 1 3allbll3) and non-minimally coupled 
models ( Pace et al.ll2014h . 


On the other hand, recent developments in observational cos¬ 
mology using high quality data including Type la supemovae (SNe 
la), cosmic microwave background (CMB), baryonic acoustic os¬ 
cillations (BAO) and Large Scale Structure (LSS), converge to 
a standard cosmological model in a spatially flat geometry with 
a cosmic dark sector usually in the form of pressureless cold 
dark matter (CDM) and dark energy (DE) with negative pres¬ 
sure, respectively, in order to interpret the observed flat rota¬ 
tion curves o f spiral galaxies and the accelerated expansion of 
the Universe jRies s et al.lll9^: IPerlmutter et al.lll999l: Jaffe et al.l 


2001J; Riesset_dj2004 ; ^bHn^ket_alj200j;|EisOTSteinet alj2005l: 

Riessetal ll2007l:lHo et al.ll2008l:IPercival et al.ll20ld:ljarosik et"^ 


201 ll : Komatsu et al.l 201 ll : Ade et alj|2013allt| ). On the basis of 
the Planck experiment results i Ade et al.ll2013ai 9). DE amounts to 
~ 68%, CDM and usual baryons to ~ 27% and 5% of the total 
energy budget of the Universe, respectively. 

Einstein’s cosmological constant A with time independent 
Equation-of-State (EoS) parameter wa = —1 is the first and 
simplest candidate to describe DE. Although A is the sim¬ 
plest model, it suffers from severe theoretical and conceptual 
probl ems: the fine-tuning and t he cosmic coincidence prob- 
lems dSahni & Starobinsk^ l200d : IWeinber3 1 19891 : ICarrollI 1200 ll : 
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[Peebles & Ratrall2003l : l^dmanabh^l2003l : ICopeland et al.ll200^ . 

Moreover since A is constant in space and time, it does not clus¬ 
ter and has a negligible contribution to the energy density budget 
of the Universe at high redshifts, so that it affects the evolution of 
structures in the Universe at z ^ 1. 

In order to solve or at least alleviate the theoretical prob¬ 
lems of the ACDM Universe, a wealth of dynamical DE 
models with a time varying equation-of-state (EoS) parame- 


ing quintessence 

Wetterichl 1 19881), phantom jCaldwelll l2002l: 

Noiiri & Odintsov 

2003allbt). k-essence (Armendariz-Picon et al.l 

200ll:IChiba et al.l 

2000l). tachvon dSenl 20021; 

Padmanabhanll 2 OO 2 I) 

and dilaton Joasnerini & Venezianol |2002|; 

Arkani-Hamed et al.l 


ical DE models which have been discussed extensively in literature. 

Other dynamical DE models which can successfully inter¬ 
pret the current accelerating Universe are constructed on the ha- 
sis of quantum gravity theories. Models such as holographic dark 
energy (HDE) models jHofava & MinidEoOd: lThomasll2Q02h and 


agegraphic dark energy (ADE) models Cail 2007tl are derived in 


the framework of quantum gravity, hy introducin g a new degree of 
freedom or by mod i fying the theory of gravity jHofava & Minid 
l200d : l^omadbool : ICafcoO^ . Recently, ghost dark energy mod¬ 
els have attracted a lot of interests in the category of dynamical 
DE models (see sectionj^for a in depth description of this class of 
models). 

Dynamical DE models not only alleviate the theoretical prob¬ 
lems affecting the cosmological constant A, but also, like pressure¬ 
less matter, possess fluctuations. Hence these models directly af¬ 
fect: 

i) the dynamics of the background cosmology through the modifi¬ 
cation of the Huhble parameter, 

ii) the matter power spectrum and large scale clustering via their 
fluctuations. 

Two parameters are important and give a significant influence to DE 
perturbations on the matter power spectrum and large scale struc¬ 
ture formation. The first is that the EoS parameter of DE should 
be different from —1. The second is the effective sound speed Ce 
which connects the density and pressure perturbations according to 
5p = c^Sp. In fact, the sound horizon of DE is given by CaH~^. 
DE can cluster only on scales larger than the sound horizon and 
smaller than the gravitational horizon, CeH' "^ < A < H-^. Here 
we consider two different cases (in units of the speed of light, c) 
usually discussed in literature: Ce = 1 (smooth DE) and Ce = 0 
(clustering DE). In the first case, the sound horizon of DE is equal 
to the Hubble length, so that DE perturbations occur on scales equal 
to or larger than the Hubble horizon and are therefore negligible on 
sub-Horizon scales. In the second scenario, the sound speed is very 
small compared to the speed of light c. Therefore the sound hori¬ 
zon of DE is significantly smaller than the Hubble horizon. In this 
case one can consider DE perturbations with a wavelength larger 
than the sound horizon and smaller than the gravitational Hubble 
horizon which can grow in a similar fashion to matter perturbations 
which growth under gravitational instability llAbram o et al.| [2009l : 

lApplebv et alj2oi^:IPace et aD2014l:iMehrabi et alj20l4). 

The SCM has been investigated and improved in clustering 
DE models in order to study how DE pertur bations impact struc¬ 
ture formation in the hi ghly non-linear regime l lAbramo et al.l200^ . 
1 20091 : IPace et alll2014l) . In this work we extend the SCM in ghost 
DE models by taking into account the perturbations of the DE fluid. 
We study the evolution of overdensities in ghost DE models and ob¬ 
tain the linear overdensity threshold for collapse Jc as well as the 


virial overdensity Avir parameters. We show how these quantities 
are affected by DE clustering in ghost DE models. We then follow 
the SCM in clustering ghost DE Universes by calculating the mass 
function and the cluster number count with the Press-Schechter for¬ 
malism. 

The paper is organized as following. In sectionl^we introduce 
the ghost DE models and describe the evolution of background 
cosmology in these models. In section the non-linear SCM in 
clustering ghost DE models is presented. In section |4] we compute 
the predicted mass function and cluster number count in ghost DE 
models using the Press-Schechter formalism by taking into account 
the clustering DE sector. Einally we conclude and summarise our 
results in section |5] 


2 BACKGROUND COSMOLOGY IN GHOST DARK 
ENERGY MODELS 

The origin of ghost fields traces back to Veneziano ghosts which 
have been proposed to find a solution to the (7(1) problem in 
the low energy effective QCD theory i Wnezian^ 19791: WittenI 
1 19791 : iKawarabavashi & Ohtalll980l : iRosenzweig et alJll98(ll) . Al- 

though in this formalism the ghost field has no contribution in 
the flat Minkowski spacetime, in curved spacetime it has a small 
energy density proportional to AqqqTT, wher e H is the H ubble 
parameter and Aqcd is the QCD mass scale ( IOhtall2oIIh . With 
Aqcd ~ 100 MeV and H ^ 10“®® eV, the quantity Aqqj^H 
gives the right order of mag nitude of the observed energy density 
of DE (~ 3 X 10“® eV)^ jOhtall201ll) . This small vacuum en¬ 
ergy density can be considered as a driver engine for the evolution 
of the Universe. Comparing with other theoretical dynamical DE 
models like HDE and ADE models in which one should introduce 
a new parameter or a new degree of freedom, the most important 
advantage of ghost DE models is that they come from the stan¬ 
dard model of particle physics and are totally embedded both in the 
standard model and in general relativity without introducing any 
new parameter or new degrees of freedom. This numerical coinci- 
denc e also shows that this mod el can so lve the fine tuning prob¬ 
lem JUrban & Zhitnits^l2009al in I2010al f9). As it was mentioned 
above, the energy density in g host DE model is proportional to the 
Hubble parameter ( IOhtal201l]) : 


Pde = aH , (1) 

where a > 0 is roughly of order Aqqu. From the observational 
point of view, ghost DE models can fit a whole set of cosmo¬ 
logical data including SNe la, BAO, CMB, Big Bang nucle osyn- 
thesis (BBN) and Hubble parameter data dCai et al] |201lh . The 
cosmological evolution of th e mode ls has been investigated in 
dUrban & Zhitnits^l2009al lhl I2010al l9l who stated that the Uni¬ 
verse begins to accelerate at redshift around z ~ 0.6. 

The Friedmann equation for a Universe containing pres¬ 
sureless dust matter and DE in a flat geometry described by a 
Friedmann-Robertson-Walker (FRW) metric is given by 

^" = ^(pm+Pde), (2) 

where pm and pde are the energy density of the pressureless dust 
matter and DE components, respectively, and H is the Hubble pa¬ 
rameter. We use equation m for the energy density of the DE com¬ 
ponent and insert it into equation ^ in order to obtain the Hubble 
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parameter in ghost DE cosmologies 


„ 47rG 
H^—a + 



(3) 


In terms of the dimensionless energy density = 

8TTGpin,o/(3HQ) and redshift parameter z = l/a — 1, the above 
Hubble equation becomes 

H{z) = Ho (^K + , (4) 

where k = (1 — f2in,o)/2. 

In the above mentioned studies, the energy density of ghost 
DE is assumed to be proportional to the Hubble parameter via equa¬ 
tion O. However, the energy density of the Veneziano ghost field 
in QC D theory is generally of the form H + 0{H^) dZhitnitsk^ 
I 2 OI ih . Although in ghost DE models only the leading term H was 
assumed, it has been shown that the sub-leading te rm can also 
be im portant in the early evolution of the Universe l lMaggiore et alj 
[2 oI 3) . A ghost DE model with sub-leading t erm i s usu- 
ally called generalised ghost dark energy model. ICai et alJ 1 I 2 OI 2 II 
showed that generalised ghost DE models result in better agreement 
with observations compared to ordinary ghost DE model s. The en¬ 
ergy d ensity of a generalised ghost DE model is given by dCai et al.l 

[ 2 ^ 

Pde = aH + , (5) 


where a and P are the constants of the model. Inserting equation © 
into equation IS, the Hubble parameter becomes 


drrG 

37 



SttG 

—-pmOCt 

37 


( 6 ) 


where 7 = 1 — 87 rG/ 3 / 3 . Setting 7 = 1 , generalised ghost DE 
model reduces to ordinary ghost DE model {P = 0), as expected. In 
terms of the dimensionless matter energy density Hm,o and redshift 
parameter z, the Hubble equation in generalised ghost DE models 
is given by 


H{z) = Ho 



_2 , f2ni,o(l + ZP 


+ 


(V) 


where k = (1 — f 2 m,o/ 7 )/ 2 . 

The conservation equations for pressureless dust matter and DE at 
the background level are 


Pm “t“ 3fTpm — 0 , 

Pde T 3fT(l A tnde)pde ~ 0 , 


(8) 

(9) 


where the dot is the derivative with respect to cosmic time and Wde 
is the DE EoS parameter. 

Taking the time derivative of Friedmann equation © and us¬ 
ing equations and © as well as the critical density pc = 
3H^/{SttG), we obtain 

^ = -^(1 + f^deWde) , (10) 

where fide is the dimensionless density parameter of the DE com¬ 
ponent. Differentiating equations 0 and 0 with respect to time 
and inserting the results in the conservation equation for DE (equa¬ 
tion © and also using equation 01 , the EoS parameter for ghost 
DE and generalised ghost DE models are 


Wde ( 2 ) 
Wde ( 2 ) 


1 

Hde(2) -2 ’ 

1 fide( 2 ) 7 

flde(2) (1 - flde(2) + 7 ) 


( 11 ) 

( 12 ) 


Setting 7 = 1 , equation GU reduces to GD as expected. We now 
calculate the equation of motion for the energy density of DE in 
ghost DE and generalised ghost DE models. Taking the time deriva¬ 
tive of fide = pde/pc and using equations 0 and dlOb and finally 
changing the time derivative to a derivative with respect to cosmic 
redshift z, we have 

C(flde(2) 3flde(2) \ / M 

= -WT7) ’ (13) 

for both ghost and generalised ghost DE cosmologies. We see that 
the evolution of the DE density in ghost DE cosmologies via equa¬ 
tion G3 depends on the EoS parameter of the models according to 
the relations dl lb and ( 112b . We solve the system of coupled equa¬ 
tions 0 , (GD, Gill and O, Gil, G1 in order to calculate the 
evolution of the Hubble and EoS parameters as well as the energy 
density of DE in ghost and generalised ghost DE cosmologies, re¬ 
spectively. To fix the cosmology, the present values of matter den¬ 
sity and DE density parameters are chosen as: fim.o = 0.27 and 
fide.o = 0.73 in a spatially flat Universe. The present Hubble pa¬ 
rameter is Ho = 70 km/s/Mpc. 

In figure 0 we show the evolution of the EoS parameter Wde 
(top panel), dimensionless Hubble parameter E = H/Hq (middle 
panel) and energy density of DE component flje (bottom panel) 
as a function of the cosmic redshift 2 for different ghost, gener¬ 
alised ghost DE and ACDM cosmological models. In the case of 
the generalised ghost DE model we choose the model parameter 
7 as 0.96 and 0.90. In all the figures of this work we indicate the 
ghost DE as GDE and the generalised ghost DE as GGDE. We see 
that the EoS parameter for ghost and generalised ghost DE models 
is always bigger than wa ~ —1 and remains in the quintessence 
regime, i.e., Wde > — 1. The Hubble parameter and the DE density 
are bigger in these models compared to the concordance ACDM 
Universe. This means that the rate of acceleration is stronger for 
(generalised) ghost models and that the relevance of the DE fluid is 
important for a longer period of the cosmic history. 


SPHERICAL COLLAPSE IN GHOST DARK ENERGY 
COSMOLOGIES 


In this section we investigate the SCM in the framework of cluster¬ 
ing ghost and generalised ghost DE cosmologies. The effects of DE 
perturbations on the evolution of matter over densities have been ex¬ 
tensively investigated in lit e rature (see, i .e. |Mota & van de Brucd 
2004 : Abramo et al. 2007 . 20081 . 20091: 0~emindli et alj 20 id : 


Basseet^ 11201 ll: iBatista & Pacell2013l: IPace et al.ll2014h. Follow 


mg 


Pace et aO ( 20141) . the fully perturbed equations for the evolu¬ 
tion of non-relativistic dust matter and dark energy 5de pertur¬ 
bations in the non-linear regime (without the contribution of shear 
and rotation) are given by 


e 


+ (1 + ^m) — — 


^de Wde^de A [1 A Wde A ^de] 


/ 2 E'\ P 3 _ , ^ r l 

0 A (-1- — ) 0 A -T-h — [flDM(5m A flde^de] 

\a E J 3a 2a 


0(14) 
0(15) 
0 (16) 


where 9 = 9/H is the dimensionless divergence of the comoving 
peculiar velocity for dust matter and DE. 
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Figure 1. Top panel: Evolution of the EoS parameter w^q. Middle panel: 
dimensionless Hubble parameter E{z). Bottom panel: DE density param¬ 
eter Qde ^ function of cosmic redshift 2 for the different cosmological 
models considered in this work. As shown in the legend, the ghost DE is 
indicated by GDE, generalised ghost DE by GGDE. The ACDM model is 
shown with the orange dashed line, the GDE model with the red short-short 
dashed line and the GGDE model with the brown (pink) short-dashed line 
for 7 = 0.96 (7 = 0.90), respectively. 


The linear evolution of overdensities at early times is 


+ ■ 


3 0 

^de U^de^de H” [l H” ti?de] 

a a 


^ + (-^- — ) 0 + ■— [riDM^m + i^de^de] 


E 


2a 


0, 

(17) 

0, 

(18) 

0 

(19) 


As in IPace et al.l j2014h . to determine the initial conditions to 
solve the above differential equations describing the evolution of 
perturbations, we find the initial value ^m,i such that at the collapse 
scale factor Uc the matter overdensity diverges, Jm —>■ 00 . The ini¬ 
tial values for the DE overdensity (5de.i ^nd peculiar velocity per- 
turbation 6i are related to (5m,i via teatista & Pacel2Q'l3l : |Pace et al.l 


|2014|): 


Ti 

<^de,i — / 0 \ '^de)<5m,i j 

(n — 3w) 

(20) 

Oi = 77'(5m,i ■ 

(21) 


In the case of an EdS model n = 1. However, in DE cosmolo- 

i ;ies it has been show n that there is a small deviation from unity 
Batista & Pacell2Ql^ . In the limiting case of non-clustering DE 
models (Ce = 1) the coupled non-linear equations (ElIIslIIU and 
system of linear equations | |17II18|[T^ . respectively, reduce to the 
following equations: 


c + 



<5' 


4 S'i 

3 1 + (5m 2a® 752 



^m(l + ^m) — 0 , 


srimo 

2 a ®£;2 


^m=0, 


( 22 ) 

(23) 


as expected (see also lPace et alj20l'^ . In the case of non-clustering 
ghost and generalised ghost DE models, we will solve equations 
( 1221123b . instead of the previous ones. 


4 RESULTS 

We first solve the system of equations ( I17t . ( 118b . jl9b for the case 
of clustering DE (Ce = 0) and ( 123b for non-clustering DE (Ce = 1) 
in order to determine the linear evolution of the overdensities 5m 
and 5de (5m in the non-clustering DE case) and obtain the linear 
growth factor D+{z) = 5m(z)/5m(z = 0) as a function of cosmic 
redshift 2 for the different cosmological models considered in this 
work. 

In figure ([3, the evolution of the growth factor normalised at 
2 = 0 and divided by the scale factor a is presented as a function 
of redshift. In the EdS model (black solid line) D+ /a is equal to 1 
at any time showing that the growth of matter perturbations 5m is 
the same at all redshifts. In a ACDM Universe (orange long-dashed 
curve) the growth factor is higher than the EdS model throughout 
its history, but falls for lower redshifts. The lowering of the growth 
factor at lower redshifts is due to the fact that at late times the cos¬ 
mological constant dominates the energy budget of the Universe 
and suppresses the amplitude of perturbations. On the other hand, a 
larger growth factor in a ACDM Universe at higher redshift shows 
that the growth of dust matter perturbations will be stronger than 
in a EdS Universe at early times. For ghost and generalised ghost 
DE models we show D+ja for two extreme cases: fully cluster¬ 
ing DE (Ce = 0) and homogeneous DE (Ce = 1) models. In the 
case of generalised ghost DE model, we choose a value of 0.96 
for the model parameter 7 . The red and brown dashed curves stand 
for clustering ghost and generalised ghost DE models, respectively. 
The blue dotted-dashed and green dotted curves show the varia¬ 
tion of D+ / a for homogeneous ghost and generalised ghost DE 
Universes, respectively. We see that the growth factor is largest 
(intermediate) for homogeneous (clustering) ghost and generalised 
ghost DE models compared to the concordance ACDM Universe. 
It is worth to mention that D+/a changes more rapidly at low 
redshifts where dark energy dominates the energy budget of the 
Universe. At early times instead, the relative change is much shal¬ 
lower. Differences are obviously still noticeable in the models. In 
particular we notice that while D+ / a for the ACDM model is al¬ 
most a constant, this is not entirely the case for the dark energy 
models analysed in this work. This is easily interpreted taking into 
account the stronger importace of dark energy at early times, as 
clearly shown in the lower panel of figure ([D- Similarly to the 
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Figure 2. Time evolution of the growth factor normalised to the present 
time and divided by the scale factor a as a function of the cosmic redshift z 
for different cosmological models. Black solid line shows the EdS model, 
the orange long-dashed line the ACDM model, the blue (green) dot-dashed 
(dotted) line the homogeneous ghost (generalised ghost) DE model, while 
the red (brown) dashed line the clustering ghost (generalised ghost) DE 
cosmology. 



ACDM Universe, in ghost and generalised ghost DE models, DE 
suppresses the growth of perturbations at low redshifts. More quan¬ 
titatively, at high redshifts, the growth of perturbations parameter 
D+ /a in non-clustering (clustering) generalised ghost DE model is 
« 18.7% (« 10.9%) larger than in the ACDM concordance model. 
These values are only midly larger than the ghost dark energy mod¬ 
els, where D+ja differs of « 14.5% (« 10%) for non-clustering 
(clustering) ghost DE models with respect to a ACDM model. 


4.1 Spherical Collapse Model parameters 

Here we evaluate the two characterising quantities of the spher¬ 
ical collapse model: the linear overdensity parameter 5c and the 
virial overdensity Avir for ghost and generalised ghost DE cos¬ 
mologies. The linear overdensity 5c together with the growth factor 
D+ is important in or der to evaluate the mass function in the Press 
Schechter formalism jPress & Schechteilll974l : ISond et aljfl^: 
ISheth & Tormenll2002n . The virial overdensity Avir is used to cal¬ 
culate the size of spherically symmetric halos with given mass M. 

We s earch the initial conditio ns 5m,i following the general ap¬ 
proach in IPace et al.l | |2010L|20I^ and then we use equations 1201 
and HD to determine 5de,i and 9i. Once the initial conditions are 
found, we then solve the system of equations GD, (DU and d in 
order to obtain 5c = 5ni(z = Zc)- In the non-clustering case, we 
solve equation 123b to calculate 5c. 

The collapse redshift Zc is defined as the redshift a t which the 
matter overdensity tends to infinity 5m —oo (see also lPace et al.l 

l2010ll20ll[20ll . 

In figure ([D the evolution of the linear overdensity 5c as a 
function of the collapse redshift Zc is presented for different mod¬ 
els: EdS, ACDM, homogeneous and clustering ghost DE and clus¬ 
tering generalised ghost DE. We refer to the caption for line styles 
and colours of each model. Analogously to the previous section, 
we use 7 = 0.96 in the case of generalised ghost DE models. In all 
DE models, the linear overdensity parameter approaches the fidu¬ 
cial value in the EdS Universe 5c « 1.686 at high enough red¬ 
shifts, as expected. In fact at high redshifts, the Universe is domi¬ 
nated by pressureless dust matter and the effects of DE on the sce¬ 
nario of structure formation is negligible. At lower redshifts, 5c de- 



Zc 

Figure 3. The variation of the linear threshold density contrast 5c as a func¬ 
tion of the collapse redshift for the different cosmological models analysed 
in this work. Line styles and colours are as in Fig. [2] 


creases and deviates from the EdS limit. The important point to 
note is that in clustering ghost and generalised ghost DE models, 
the behaviour of 5c is more similar to the ACDM Universe with 
respect to the homo geneous models. This result is in agreemen t 
with what found by Pace et al. I Hold) an dlSatist a & Pacel HolSh 


for early dark energy models and IPace et al.l Hoidh for clustering 
quintessence and phantom DE models. Quantitatively we see that 
the differences between homogeneous ghost and generalised ghost 
DE models with the ACDM model at the present time are roughly 
4.5% and 21.5%, respectively, while in the case of clustering ghost 
and generalised ghost models the differences are very small, of the 
order of « 0.25% for both cases. As for the case of perturbations 
in early dark energy models, we can explain the fact that perturba¬ 
tions in the dark energy fluid make the model more similar to the 
ACDM model by taking into account that the Poisson equation is 
now modified and the gravitational potential is sourced also by DE 
perturbations. 

In the framework of the SCM, the virialization pro¬ 
cess of pressureless dust matter and the size of forming 


199ll: Wang & Steinhardt 19981: Mota & van de Bruck 2004 

Horellou & Bergell2005l; 

Wang & TegmarkI 

20051) and also its per- 

turbations dAbramo et al 

200l 20081. 2009 

: Batista & Pacell2013l: 


Pace et. 


(^{xjy)^, where is the overdensity at the turn-around epoch, x 
is the scale factor normalised to the tum-around scale factor and 
y is th e ratio between the viriali zation radius and the turn-around 
radius dWang & Steinhar^ll998l) . In an EdS cosmology, it is sim¬ 
ple to show that y = 1/2, « 5.6 and Ayjr ~ 178 i ndepen- 

dently of the cosmic redshift (see also lNaderi et al.ir2015h . In DE 
cosmologies, Avir depends on the evolution of the DE fluid and it 
is a redshift dependent quantity. In particular, according to whether 
DE takes part or not into the virialization process, the quantity y 
may be larger or smaller than 1/2 and the parameter Ayjr can be 
affect ed by the clustering of DE dMaor & Lahavll200^ : IPace et all 
|2014 . 


In the line of these studies, we calculate the turn-around and 
virial overdensities C, and Avir in ghost and generalised ghost DE 
cosmologies. We also determine how the clustering of the DE com¬ 
ponent can change the variation of </ and Avir in these models. Our 
results are presented in figure (|4). In the top panel, the variation 
of the overdensity at the tum-around redshift (/ is shown for the 
different DE models investigated in this work. In the limiting case 
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Zc 



Zc 


Figure 4. The variation of turn around overdensity ^ (top panel) and virial 
overdensity Avir (bottom panel) with respect to the collapse redshift Zc for 
vai'ious models considered in this work. Line styles and colours are as in 

Fig.El 


of the EdS model, C, = 5.6 independently of the cosmic time. At 
early times, tends to the critical value ^ = 5.6 representing the 
early matter-dominated era. The value of C, is larger for both clus¬ 
tering and non-clustering versions of ghost and generalised ghost 
DE models, compared to the concordance ACDM model. Differ¬ 
ences between the dark energy models and the ACDM model cover 
a relatively huge spectrum of values, according to whether we take 
into account the perturbations of the dark energy fluid. In particular 
for homogeneous (generalised) ghost DE models, differences are 
of the order of (67%) 17% while when perturbations are taken into 
account, differences drop down to (9%) 10%. This shows how the 
generalised ghost model, is much more sensitive to the inclusion of 
DE perturbations. Once again, this is due the higher amount of DE 
density at earlier time for these kind of models. Hence we conclude 
that in the context of ghost DE scenarios the perturbed spherical 
perturbations detach from the background Hubble flow with higher 
overdensities compared to the EdS and ACDM Universes. We also 
notice that the results for clustering ghost and generalised ghost DE 
models are closer to what is found for the reference ACDM model 
compared to the non-clustering case. 

In the bottom panel, the results for the virial overdensity Avir 
are presented. In the case of ghost DE models, the results are closer 
to the reference ACDM model. In these models, the difference be¬ 
tween homogeneous and clustering DE is very small. However, in 
the case of generalised ghost cosmology the differences are more 
pronounced. Differences between ghost DE models and ACDM 
model at the present time are of the order of 60% for non-clustering 
generalised ghost DE model and 7% for homogenous ghost DE 


models, roughly as for the overdensity at the turn-around radius. 
This is not surprising, since this quantity is the main ingredient 
used to evaluate the virial overdensity. When perturbations in dark 
energy are taken into account differences are only of the order of 
1.5% fo r both classes analysed. Thes e results are simi lar to what 
found in idel Ponolo et al] ( 1201 3d) and IPace et alj ( l2014h . 


4.2 Mass function and halo number density 


We know that galaxies and cluster of galaxies are embedded in 
the extended halos o f cold dark matter (CDM ). In the Press & 
Schechter formalism dPress & Schechte3ll974h . the abundance of 
CDM halos can be described as a function of their mass and a 
Gaussian distribution function expresses the fraction of the volume 
of the Universe which collapses into an object of mass M at a cer¬ 
tain redshift 2 . In this formalism, the comoving number density of 
virialised structures with masses in the range of M and M + dM 
at redshift 2 is given by 

dn{M,z) p^Q dlna{M,z) 

-^M- = ~ir dM ’ ^24) 

where a is the r.m.s. of the mass fluctuation in spheres of mass M 
and /(cr) is the mass function. The stan dard mass function in the 
Press & Schechter formalism is given by dPress & Schechteill974h 


/(<^) 



\ s!iz) 1 

/ TV a{M, z] 


2(t2(M, z) 


(25) 


Although the standard mass function provides a good general 
representation of the predicted number density of CDM halos, it 
deviates from simulations due to an over-prediction of low-mass 
objec ts and an under-prediction of h igh-mass objects at the present 
time dSheth & TormenI [19991 l2002h . Here we u se another popu¬ 
lar fitting fo rmula proposed by Sheth & Tormen dSheth & TormenI 
I1999L|20 o 3) . the so-called ST mass function: 


/sT(cr) = A 



Sc{z) 
cr{M, z) 


exp 


aSl{z) 
20-2 (M, 2) 
(26) 


where the numerical parameters are: A = 0.3222, a — 0.707 and 
p = 0.3. Putting A = 1/2, a = 1 and p = 0, the well known 
Press-Schechter mass function can be recovered as expected. In a 
Gaussian density field, the mass variance is given by 

k^P(k)W'^{kR)dk , (27) 

where R — (iM/A -k is the radius of the overdense 
spherical path at the present time, W{kR) = 3[sin(fci?) — 
kRcos{kR)\/{kR)^ is the Fourier transform of a spherical top- 
hat profile window function with radiu s R and P(k) is the linear 
power spectrum of density fluctuations l |Peebleslll993h . 

The number density of objects above a given mass at a certain 
fixed redshift is 

"<> ">=r ™ 

In this section, having at hands all the necessary ingredients, 
we can compute the predicted number density of virialised ob¬ 
jects in the Press-Schechter formalism for ghost and generalised 
ghost DE cosmologies. As first step, we use equation ( 128b for non¬ 
clustering ghost and generalised ghost DE models. In next section 
we investigate the predicted number densi ty of virialised clusters o n 
the basis of the formulation presented in Icreminelli et all ( l2010h : 
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Figure 5. Ratio of the number of objects above a given mass M for halos at z = 0 (top left), z = 0.5 (top right), z = 1.0 (bottom left) and z = 2.0 (bottom 
right) between the ghost and generalised ghost DE models and the concordance ACDM model. The blue dotted-dashed curve represents the ghost DE model. 
The purple dashed, green dotted and pink dashed curves stand for generalised ghost DE model with model parameter 7 = 0.98, 7 = 0.96 and 7 = 0.90, 
respectively. 


iBasse et ^ 1 I 2 OI ih i lSatista & Pacd ll201.3h where the total mass of 
the halos is affected by DE perturbations. 

We use the /st mass function given in equation d26b and 
choose four different redshifts: z — 0, z = 0.5, z = 1.0 and 
z = 2.0. In order to calculate a , w e follow the formula tions pre¬ 
sented in lAbramo et alj ( l2007h and iNaderi et"^ ( l2015h . We also 
adopt the ACDM model as reference model with normalization of 
the matter power spectrum aa = 0.776, in agreement with recent 
observations (??). In figure IS, we show the ratio of the number 
of objects above a given mass M between the ghost and gener¬ 
alised ghost DE models and the concordance ACDM model, in the 
case of non-clustering DE cosmologies. In the case of generalised 
ghost DE model, we choose three different values for the model 
parameter; 7 = 0.90, 7 = 0.96 and 7 = 0.98. As was discussed 
in section IS, for 7 = 1 . 0 , the generalised ghost model reduces 
to the ghost DE model. In the upper left (right) panel, the results 
are presented for redshift z = 0 (z = 0.5). The lower left (right) 
panel is for halos at redshift z = 1 (z = 2 ), as indicated in the 
legends. Due to the identical normalization of the matter power 
spectrum, we see that at z = 0 all models have the same num¬ 
ber of objects. At z = 0.5, the ghost (7 = 1) and generalised 
ghost with 7 = 0.96 and 7 = 0.98 are still giving the same num¬ 
ber of object compared to the concordance ACDM model, with 
very small differences for the large mass tail of the distribution 
{10^^Mq/H). However, for smaller values of 7 the differences be¬ 
tween the generalised ghost and the ACDM model are consider¬ 
able. For 7 = 0.90, the predicted number density for virialised 
halos at high mass tail (10^^Mq/ h) is roughly 40% lower than the 
concordance ACDM model (see top-right panel of figure (Ell. At 
higher redshifts, z = 1 and z = 2, the predicted number of struc¬ 


tures exceeds what predicted by the ACDM model, in the cases of 
ghost and generalised ghost DE models with model parameter 7 
close to unity, i.e., 7 = 0.96, 0.98. The generalised ghost model 
with low value 7 = 0.90 shows a decrement in the number of 
objects compared to other cases. We see that a major difference 
between ghost DE models and ACDM model takes place at high 
mass while all models are roughly producing the same number of 
objects at the low mass tail, as expected. From a more quantita¬ 
tive point of view, at redshift z = 1, the ghost DE model shows 
an increment in the number of structures of « 47%, while for the 
generalised ghost DE model with 7 = 0.98 (7 = 0.96) the incre¬ 
ment is « 38% (« 29%), when they are compared to the standard 
ACDM model. Surprisingly, the case with 7 = 0.90 shows a sub¬ 
stantial lack of high mass objects, of the order of « 66 % for objects 
of mass « lO^®AT 0 //i. 


4.3 Corrected mass function in clustering ghost DE models 


It is well known that in EdS cosmology y = i?vir/77ta = 1/2. 
In the EdS model, the virial overdensity can be calculated ana¬ 
lytically and the calculations lead to Avir — 178 independently 
of the redshift z. It is also well known that this value strongly 
depends on the particular background cosmology and changes in 
the presence of DE. It was shown that the virialization process of 
dark matter overdensities in the non-linear regime depends on the 
properties of DE mod els iL ahav et a illl99ll ; iMaor & Lahavll200^ ; 
ICreminelli et alJl 20 ld ; lBasse et al. 2011 ). Moreover, in clustering 
DE models, we should take into account the contribution of DE 
perturbations to the total mass of the halos fC reminel li et akll^ld ; 

ISasse et ^l201ll ; iBatista & Pacd[20l^ ; IPace et al.ll2014t) . Depend- 
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Figure 6. The ratio of DE mass to dark matter mass by using the definition 
of e via equation (ID in the top panel and equation (ID in the bottom panel. 
The red dotted-dashed curve represents the clustering ghost DE model. The 
purple long dashed, brown dashed and pink long-dashed curves stand for 
clustering generalised ghost DE model with model parameter 7 = 0.98, 
7 = 0.96 and 7 = 0.90, respectively. The orange long dashed one indi¬ 
cates the ACDM model. 


ing on the quintessence or phantom DE EoS parameter, Wde (z), DE 
can add or subtract mass to the total mass of the halo, respectively. 
The fraction of DE mass to be taken into account with respect to the 
mass of the dark matter is given by the quantity e{z) = Mde/Mm- 
The mass of dark matter Mm is defined as 

r^vir 

Mn, = / dR R\i + S^) . (29) 

Jo 

where pm is the mean density of dark matter. The mass of DE in 
the virialization process depends on what it is considered to be the 
mass of the DE component. If we only assume the contribution of 
DE perturbation, then the mass of the DE component is given by 

/■Rvir 

Ml, = 47rpde / dRR^5d,{l + 3Ce) , (30) 

Jo 

where pde is the mean density of dark energy. On the other hand, if 
we assume also the contribution at the background level (in analogy 
to dark matter), the total mass of DE is given by 

/■Rvir 

^de = 47rpde / dRR^ [(1 + SlUde) + (5de(l + Sc^)] . (31) 

Here we compute the quantity e{z) in the case of full clus¬ 
tering DE scenario. It should be noted that due to the background 
contribution in equation OH . even in the case of homogeneous DE 
scenario s, DE can add or subtract mass to the dark matter halos 
(see also lSatista & Pacell2013l : IPace et al.ll20l4) . Since we work in 


the framework of the top-hat spherical profile, the quantities inside 
the collapsing sphere evolve only in time without any spatial de¬ 
pendency. Hence on the basis of the definition in equation OOt we 
have 


e{z) 


Hde(2^) ^de 

1 + <5m ’ 


and under the same assumption for equation (ED 


(32) 


/ Dde(^) 1 3mde(^) “t" ^de /oi, 

<^( 2 ) = w / ( -;—- • (33) 

Din(. 2 ) 1 -f <5m 

In figure (D we show the evolution of e(z) according to equa¬ 
tions d32l l (top panel) and l l33b (bottom panel) for both ghost and 
generalised ghost DE models. We refer the reader to the caption for 
line styles and colours. We notice that according to equation d32b . 
both ghost and generalised ghost DE models give a positive con¬ 
tribution to the total mass of the halos (e{z) > 0). The quantity e 
grows with redshift z meaning that the contribution of DE on the 
total mass is larger for halos which virialize at lower redshifts. This 
is explained by taking into account the fact that at late times dark 
energy perturbations arise, while being negligible at earlier times. 
In particular, the generalised ghost DE model gives a higher con¬ 
tribution of DE to the total mass of halos compared to ghost DE 
model. We also see that a lower value of the model parameter 7 
results in a higher value of e(z). On the other hand, using the def¬ 
inition in equation d33b , the behaviour of e is different. Looking at 
equation ED, we see that the sign of the quantity 1 + 3mde + (ide 
determines how DE contributes to the total mass of the halo. During 
the history of the halo, when 3u;de + ^de < — T the contribution of 
the DE mass is negative and therefore it lowers the total mass of the 
halos. The other important point to note is that according to the def¬ 
inition of e via equation d33b . one can consider the DE mass even 
for homogeneous DE models (Jde = 0). Hence in the bottom panel 
of figure ( 6 ) we show £ for the ACDM model (rcde = — !)• We also 
see that for the ACDM model e is always negative. On the other 
hand, when 3u;de + 5de > — 1, the contribution of DE is positive. 
This is the behaviour taking place for ghost and generalised ghost 
DE models. Comparing the top and bottom panels in figure l|D, 
we see that for ghost and generalised DE models the values of e 
according to the second definition are lower than equation l l32b . 

We now compute the number density of virialized halos in the 
pr esence of DE mass co rrecti on. Follow i ng the procedure outlined 
in lSatista & Pace! ( l2013h and IPace et alj ( l2014h . in the presence of 
DE perturbations, we change the mass of halos as M —>■ M(1 — e). 
In this case the correc ted mass function can be redefined as follows 
(Batista & Pacel2013tl 


dnc{z,M) pmo dlna{M,z) , , 

dM =-M(1-£) dM 

where /(ct) is given by equation l l26b . It should be noted that the 
clustering of the DE component can also change the mass function 
/(a) by changing the quantities S, and (j(M, z). In figure ([7) we 
show the relative number density by using the corrected mass func¬ 
tion formula in equation d34b and the definition of e in equation 
ED- Analogously to figure 0 . we choose four different redshifts 
z = 0, 0.5,1 and z = 2. We see that due to the contribution of 
the DE mass, the relative number density between ghost and gener¬ 
alised DE models and ACDM model ric/riACDM is equal or larger 
than one. In the top panels, for 2 = 0 and 2 = 0.5, all ghost, 
generalised ghost and ACDM models generate the same number 
of objects at the low-mass end. However, at high masses, gener¬ 
alised ghost DE models with smaller value of the model parameter 
7 produce more objects. At higher redshifts, 2 = 1 and 2 = 2, 
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Figure 7. Ratio of the number of objects above a given mass M for halos at 2 = 0 (top left panel), 2 = 0.5 (top right panel), 2 = 1.0 (bottom left panel) and 
2 = 2.0 (bottom right panel) between the ghost and generalised ghost DE models and the concordance ACDM model by using the mass definition in equation 
(32). The red dotted-dashed curve represents the clustering ghost DE model. The puiple long dashed, brown dashed and pink long-dashed curves stand for 
clustering generalised ghost DE model with model parameter 7 = 0.98, 7 = 0.96 and 7 = 0.90, respectively. 


one can see that ghost and generalised ghost DE models have more 
objects for both the low mass and the high mass tail compared to 
the concordance ACDM model. Comparing figures O and 0> we 
conclude that for all the redshifts here considered, the predicted 
relative number density of halos in the presence of DE corrections 
to the halo mass is higher than what was found for homogeneous 
DE scenarios. In particular, the differences are larger for high-mass 
objects. This feature of ghost a nd generalised ghost D E models is 
in agreement with the results of iBatista & Pace! ( l2013h for inhomo¬ 
geneous Early Dark Energy models. To compare the models quan¬ 
titatively, we restrict our analysis to 2 = 0. While non-clustering 
DE models and the ACDM model have the same number of ob¬ 
jects (top-left panel of figure ((3), clustering DE models deviate 
from the ACDM model particulary at the high-mass tail (top-left 
panel of figure (7)). We see that clustering generalised ghost DE 
model with 7 = 0.90 doubles the values obtained for the ACDM 
model, while for 7 = 0.96 the increment in the number of objects 
is « 61%. Finally, for 7 = 0.98, differences are of the order of 
50%. As expected, the increase in 7 corresponds to a decrease in 
the differences with respect to the ACDM model at the high mass 
tail (of the order of 40% for the ghost dark energy model). 


5 CONCLUSIONS 

In this work we studied the non-linear evolution of structure for¬ 
mation in ghost and generalised ghost dark energy models within 
the framework of the spherical collapse model. Dark energy mod¬ 
els are described by a constant or a time varying equation of state 
w{a) whose functional form is usually purely phenomenological. 


The advantage of dynamical dark energy models with respect to 
the ACDM model is that to alleviate its theoretical problems. An¬ 
other interesting aspect is the presence of fluctuations in the dark 
energy fluid that can have an important and characteristic impact on 
non-linear structure formation. Ghost dark energy models are the¬ 
oretically well motivated and trace their origin back to the studies 
of the low energy effective PCD theory dVeneziandl 19791: WittenI 

' itii ' 


Il979l : l^warabavashi & Ohtal98(]l : lRosenzweig et aijl98C ). Here 

we concentrate on the cosmological implications of such models 
and present results that can be useful from an observational point 
of view. 


We start our analysis by studying the effects of the modifica¬ 
tion of the background expansion history of the growth factor and 
on a second step we take into account also the perturbations in the 
dark energy component. We showed that the linear growth factor 
is sensitive to the details of the model considered. In particular, 
when dark energy is homogeneous, differences with respect to the 
ACDM model are very pronounced and can be up to 50% at rel¬ 
atively high redshifts. Interestingly, when dark energy fluctuations 
are taken into account, differences become smaller and the mod¬ 
els studied get closer to the ACDM model. This is exactly what 
happens for early dark energy model and it applies to this class 
of models as well, being the amount of dark energy at early times 
significantly more important (see bottom panel of figure (Di- 

Being the equations governing the evolution of the linearly 
extrapolated overdensity parameter 5^ identical to those of the lin¬ 
ear growth factor, results are similar to what found for the linear 
growth factor. Once again when the perturbations in the dark en¬ 
ergy fluid are included in the calculations, the models examined 
resemble closer the ACDM model. This is easily explained by tak- 
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ing into account that when dark energy perturbations play a role, 
Poisson equation is modi fied and this makes th e model more simi¬ 
lar to the ACDM one Csee lBatista & Pacel2013L for a more detailed 
discussion). Similar conclusions can be drawn for the non-linear 
virial overdensity Ay. 

While the growth factor can be inferred from the determina¬ 
tion of the matter power spectrum on linear scales at a given redshift 
or by measuring the quantity fas{z), where / = d\n D+/dlna 
and as{z) = D+{z)as{z = 0), the linear extrapolated overdensity 
5c is a purely theoretical quantity and it would be very hard to infer 
it from observations. Nevertheless they can be combined together 
and be the basic ingredients to evaluate the mass function and the 
number of halos above a given mass. When keeping the dark en¬ 
ergy fluid smooth and important only at the background level, we 
showed that both the ghost and generalised ghost models give an 
excess of objects with respect to the ACDM model, except for the 
generalised ghost model with lower value 7 = 0.9. As expected, 
differences are very pronounced only for the high mass tail of the 
mass function, while at low masses the models are practically in¬ 
distinguishable. 

One might consider that if dark energy can clump, then the 
total mass of the halos should be affected. By defining the contri¬ 
bution of dark energy limited to that of the clumps, we showed in 
the top panel of figure nil that the dark energy mass can be at least 
10 % of the dark matter mass at z = 0, with differences growing 
with decreasing the parameter 7 and decreasing at high redshift 
where the contribution of dark energy is less important. Note that if 
the amount of dark energy at early times is significantly bigger than 
in the ACDM model, then the dark energy mass can still be an ap¬ 
preciable fraction of the dark matter mass. By taking into account 
this correction and evaluating the number of objects above a given 
mass, we showed that the dark energy models we studied in this 
work predict much more objects with respect to the ACDM model 
at all redshifts. This could be a good test to compare theoretical 
predictions with observations. 
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